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INTRODUCTION 


Because  growth  of  very  small  tumors  has  not  been  studied  and  documented  routinely  in  detail, 
this  project  has  increasingly  focused  on  the  detection  of  micro  metastases  containing  fewer  than  10^ 
cells.  The  detection  of  such  small  tumors  would  constitute  a  major  advance  in  diagnosis  if  extended 
to  the  clinical  environment. 

In  our  progress  to  date  we  have  demonstrated  our  capability  to  detect  sub-millimeter  metastases 
in  mice  by  light  emission  using  Light  Emission  Tomography  (LET),  a  technique  based  on 
bioluminescence  of  cancer  cells  transfected  with  luciferase,  and  the  ability  to  detect  metastases  in 
the  entire  body  of  a  mouse.  We  have  characterized  the  procedures  which  in  our  hands  result  most 
consistently  in  metastases  at  different  sites. 

We  have  begun  assessment  of  perfusion  using  fluorescence  imaging  and  will  add  Nuclear 
Medicine  methods.  The  question  we  are  currently  investigating  is,  can  we  reliably  detect  small 
tumor  foci  using  techniques  that  can  be  more  readily  translated  to  the  clinical  environment?  Our 
technical  goal  is  to  compare  detection  of  tumor  cells  to  signals  obtainable  by  Nuclear  Medicine 
techniques. 

Our  technological  advances  in  the  simultaneous  use  of  Single-photon  Emission  Computed 
Tomography  (SPECT)  have  resulted  in  a  new  micro-SPECT  based  on  cooled,  electron-multiplied 
Charge-Coupled  Devices  (EMCCDs)  with  which  we  are  starting  imaging  experiments.  Ultimately 
we  plan  to  use  bioluminescence  together  with  SPECT  to  detect  areas  where  tumor  cells  proliferate 
and  establish  metastases  and  to  identify  the  relative  times  at  which  the  images  first  appear. 
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BODY 


Light  Emission  Tomography. 


Bioluminescence  imaging  is  fast  and  particularly  useful  as  a  screening  technique,  validated  as  a 
high-throughput,  quantitative  modality  for  assessing  tumor  burden,  investigating  in  particular  the 
consequences  of  different  methods  of  administration  of  luciferin  *.  The  simplest  form  of  imaging  is 
non- volumetric,  but  has  several  shortcomings  its  ineffectiveness  in  precisely  localizing  cell 
populations  hampers  its  use  as  a  visualization  tool,  while  uncertainty  in  estimation  of  tumor  burden 
at  deep  sites  impairs  its  use  as  a  quantitative  tool..  To  obviate  these  shortcomings  while  retaining 
the  temporal  advantages  of  planar  imaging,  we  have  built  a  compact  multi-camera  system  to 
efficiently  perform  volumetric  imaging.  We  have  chosen  back-illuminated  CCDs  (originally 
Scientific  Imaging  Technologies  Inc.,  Tigard,  OR,  now  replaced  with  equivalent  CCDs  by  e2v 
technologies  inc.,  Elmsford,  NY)  because  of  their  high  quantum  efficiency  over  a  broad  range  of 
wavelengths  (50-95%  between  390  and  930  nm). 

Currently  the  system  has  four  cameras,  but  the  number  is  limited  only  by  cost,  complexity  and 
time  considerations  and  can  be  increased.  The  four  cameras  (Figure  1)  record  simultaneously 
multiple  views  of  the  animal  in  order  to  assess  the  spatial  distribution  of  the  emitted  light. 


Four-camera  system  for  3D  imaging.  At  each 
position  the  cameras  record  first  an  image  of  the 
animal’s  surface  using  reflected  light  and  then  a  BLI 
image.  The  two  sets  of  images  are  combined  for 
reconstruction.  The  cameras  are  then  rotated  in  18- 
degree  steps  using  the  gantry  shown  in  central  panel  to 
obtain  sufficient  data  for  high-resolution  3D  imaging. 
In  five  steps,  each  camera  reaches  the  position 
previously  occupied  by  an  adjacent  camera  in  order  to 
follow  light  intensity  changes. 


A  computer-controlled  rotating  gantry  allows 
imaging  at  intermediate  angular  position  to  give  the 
extended  angular  coverage  necessary  for  3D  reconstruction  in  combination  with  means  we  have 
devised  to  account  for  its  time  variation.  For  each  camera  and  angle  of  rotation,  we  obtain  both  an 
image  of  the  animal  using  external  illumination,  and  an  image  of  internally  emitted  light.  The  image 
of  the  animal  is  reconstructed  digitally  in  the  same  system  of  coordinates  as  BLI  and  is  crucial  to 
the  process  of  3D  image  reconstruction  in  our  non-contact  imaging  system.  All  cameras  in  the 
system  are  calibrated  so  that  the  resulting  images  can  be  digitally  combined.  Importantly,  the 
captured  light  intensity  changes  by  an  order  of  magnitude  depending  on  the  angle  of  detection  even 
though  the  source  intensity  remained  constant. 

Our  device  overcomes  the  problems  of  excessive  imaging  time  with  a  single  camera  and  the 
confounding  fact  that  light  emission  does  not  remain  constant,  but  changes  in  time  following 
administration  of  luciferin  {luc).  The  kinetics  of  light  emission  depend  on  numerous  factors,  such  as 
the  location  of  the  cells,  physiology,  the  method  of  luc  injection  {e.g.,  IP,  IV,  SQ  or  IT),  and  the 
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type  and  concentration  of  anesthetic  With  the  four  cameras  in  quadrature,  proceeding  in  18° 
steps,  each  camera  repeats  the  measurement  of  an  adjacent  camera  after  five  rotations  (6  images)  to 
detect  changes  in  light  emission. 

A  multi-camera  imaging  system  is  unique  in  that  with  sufficient  cameras  it  can  capture  the  light 
emission  from  the  entire  animal  surface  simultaneously.  In  addition  to  reducing  the  data  acquisition 
time,  this  innovative  approach  gives  an  imaging  system  that  is  more  robust  in  handling  luciferin 
kinetics  and  other  time-dependent  phenomena  that  may  be  present  in  the  imaging  data. 

The  standard  method  for  correcting  tomographic  data  using  rotating  camera  geometry  has  been 
to  scale  surface  intensities  obtained  at  different  time  points  based  on  an  interpolation  curve 
determined  from  images  taken  from  the  same  angular  direction  at  the  beginning  and  end  of  the 
experiment.  The  effect  is  depicted  in  the  figure  below,  which  shows  measured  and  corrected  data 
for  a  single  view. 


Time  PI  [min] 


Time  course  of  light  emission  after  injection  of  a  bolus  of  luciferin.  After  reaching  a  peak  at  the 
site  of  the  tumor  in  10  min  there  is  a  continuous  decay  of  intensity  which  falls  to  ‘A  of  its  peak  value 
over  10  minutes. 

Inherent  to  this  method  is  an  assumption  that  the  luciferin  kinetics  for  all  bioluminescent  light 
sources  in  the  animal  body  is  identical.  However,  we  know  from  pharmacology  and  drug  delivery 
that  luciferin  delivery  will  in  fact  be  heterogeneous  and  occur  at  a  tissue-dependent  rate. 

To  illustrate  the  possible  effects  of  site-dependent  light  emission  dynamics,  we  performed  an 
experiment  using  2  point  sources  located  at  different  depths  in  a  cylindrical,  turbid  phantom,  with 
one  source  decaying  twice  as  fast  as  the  other  (see  below). 
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Surface  light  intensity  for  2  point  sources  with  different  dynamics  in  a  turbid  phantom,  (ieft)  projection  obtained 
using  single  rotating  camera  and  corrected  using  interpolation  between  baseline  and  finai  images;  (middie) 
actuai  projection  collected  using  synchronous,  multi-camera  optical  imaging  system;  (right)  difference  between 
images,  showing  that  a  single  camera  system  is  susceptibie  to  errors  in  data  coilection  for  heterogeneous 
luciferin  kinetics. 

The  image  from  the  single  rotating  camera  (left  panel)  has  been  corrected  for  light  emission 
dynamics  using  the  interpolation  method  above.  In  the  middle  panel,  the  true  surface  intensity,  as 
captured  by  a  multi-camera  system  with  complete  surface  coverage,  is  shown.  The  rightmost  panel 
illustrates  the  error  in  data  acquisition  inherent  to  the  single  camera/rotating  view  approach,  shown 
by  subtracting  the  images  obtained  by  the  two  methods.  In  this  case,  the  single  camera  interpolation 
underestimates  the  surface  intensity  from  the  deeper  source.  This  would  undoubtedly  lead  to  errors 
in  depth  and  intensity  in  the  image  reconstruction.] 

In  summary,  we  have  tested  light  emission  from  various  sites  within  the  body  to  determine  the 
circumstances  under  which  we  could  reliably  detect  metastases.  We  developed  means  to  calibrate 
light  intensity  in  terms  of  photons/unit  area/sec  and  to  evaluate  the  light  intensity  emitted  by  the 
different  cell  lines  used  in  our  studies.  For  deeper  metastases  the  light  intensity  is  attenuated  and 
imaging  times  must  be  consequently  extended. 

Changes  of  emission  intensity  with  time  after  injection  of  luciferin.  The  same  head  interrogates 
the  same  field  of  view  at  regulated  intervals  (squares  in  figure)  and  thus  our  multiple  head 
configuration  allows  the  continuous  assessment  of  changes  in  light  emission.  As  the  different 
cameras  have  adjoining  fields  of  view  we  can  monitor  these  changes  and  correct  their  effect  on 
images. 

With  this  technology  we  were  able  to  detect  in  vitro  with  the  same  instrument  extended 
distributions  with  500  cells  and  foci  estimated  at  10-30  cells  when  the  cells  had  sufficient  emission 
intensities.  The  establishment  of  reliable  models  of  metastasis  to  different  sites  is  a  prerequisite  for 
our  research.  To  select  models  which  work  repeatably  in  our  hands,  we  performed  experiments  with 
different  implantation  models. 
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KEY  RESEARCH  ACCOMPLISHMENTS 


Comparison  of  Cancer  Cell  Implantion  Methods 

Bioluminescence  imaging  experiments  were  performed  using  MDA-MB  23 1  Luc  breast  cancer 
tumor  cells  using  different  implantation  techniques:  orthotopic,  IV,  and  intracardiac.  The  success 
rate,  the  time  required  to  form  metastatic  tumor,  and  the  location  of  metastases  were  measured. 
Bone  is  a  common  site  of  metastasis  in  human  cancer  and  an  animal  model  that  produces  such 
metastases  is  desirable  for  research  and  central  to  the  theme  of  our  project. 

1 -Orthotopic  tumor  Implantation 

In  the  case  of  the  orthotopic  implantation  40%  of  the  animals  developed  primary  tumors  and  none 
developed  metastases  while  the  primary  tumor  progressed.  No  metastases  were  detected  in  a  1 5 
days  interval  after  primary  tumor  excision  in  the  primary  areas  of  interest,  lungs  and  bone. 


Figure  1.  MDA-MB-231  LUC  orthotopic  tumor.  No  metastases  were  detected  by  15  days  after 

primary  tumor  removal. 


2-IV 

The  IV  inoculation  produced  metastases  in  the  lungs  of  the  animal  in  60%  of  the  cases.  On  average 
the  lung  tumors  could  be  detected  12  days  post  injection  and  they  grew  at  a  relatively  slow  rate. 


Figure  2.  Lung  metastases  developing  after  IV  injection  with  10^  stably  luciferase  transfected 
MDA-MB-23 1-LUC.  They  could  be  monitored  and  quantified  by  LET  for  an  average  of  90  days. 

The  rest  of  the  animals  (40%)  developed  metastases  in  the  long  bones  of  the  hind  legs  and  in  the 
abdomen.  Bone  is  a  common  site  of  metastasis  in  human  cancer,  thus  this  specific  animal  model  is 
very  valuable  for  research  (Fig.  3). 
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Day  12 


Day  33 


Day  47 


Figure  3.  Femoral  metastases  following  IV  injection  with  10  stably  luciferase  transfected  MDA- 


MB-231-LUC. 


3 -Intracardiac 

A  cohort  of  8  mice  were  implanted  with  MDA-MB-23 1-LUC  cells  by  injecting  them  in  the  left 
cardiac  ventricle.  All  the  animals  developed  metastases,  most  of  them  being  localized  in  the  lungs 
(Fig.  4).  The  growth  rate  of  the  lung  tumors  was  consistent  with  that  of  the  IV  implantation  model. 


$ 
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Day  12 


Day  14 


Day  33 


Figure  4.  Lung  metastases  developing  after  intracardiac  injection  with  10  cells 

One  animal  (1/8)  developed  metastases  in  the  femur  and  in  the  abdominal  cavity  (Fig.  5). 


Figure  5.  Long  bone  and  abdominal  metastases  developing  after  intracardiac  injection  with  10^ 

stably  luciferase  transfected  MDA-MB-23 1-LUC. 

It  is  important  to  validate  the  LET  method  against  other  established  methods. 
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Human  breast  tumor  cells  (10*’  stably  luciferase  transfected  MDA-MB-231-luc)  were  injected 
IV  in  a  nude  mouse.  The  animal  was  imaged  repeatedly  by  LET  and  MRI  over  several  weeks  as  the 
lung-colonizing  metastases  appeared  and  progressed  in  size  and  number.  Following  D-luciferin 
injection  (450  mg/kg,  SQ)  in  the  anesthetized  mouse,  a  set  of  20  images,  18°  apart,  was  obtained 
using  1  min  exposure  starting  5  min  post-injection.  Data  reconstruction  provided  a  3-D  model  of  the 
lung  tumors.  The  algorithm  was  also  applied  to  images  obtained  using  light  reflected  from  the  skin 
under  external  illumination  and  the  two  data  sets  combined.  The  imaging  time  varied  from  300 
s/angular  position  for  the  first  imaging  sessions  that  were  characterized  by  low  bioluminescent 
signal,  to  30  s/angular  position  for  the  later  sessions  when  the  bioluminescent  signal  was  large.  MRI 
scans  covering  the  chest  of  the  mouse  were  acquired  on  a  4.7  T  Varian  scanner  using  a  home  built 
respiratory  gating  unit.  We  obtained  contiguous  proton  density  weighted  spin-echo  MR  coronal 
slices  with  the  following  parameters:  TE  =  12  ms,  FOV=  3.2  cm  X  6.4  cm,  slice  thickness  =  1  mm, 
matrix=  64  X  128  zero  filled  to  128  X  256,  4  averages.  The  LET  images  first  detected  a 
bioluminescent  signal  22  days  after  cell  implantation,  and  revealed  the  growth  and  spread  of  the 
lung  metastases  at  weekly  intervals.  Tumors  in  the  lungs  were  first  detected  by  MRI  46  days  after 
inoculation.  The  smaller  tumors  from  the  left  lung  shown  on  the  LET  images  were  still  not  observed 
by  MRI.  The  position  and  relative  sizes  of  the  tumors  are  consistent  between  the  two  imaging 
modalities.  Nevertheless,  LET  was  able  to  detect  the  tumors  17  days  earlier  later  than  MRI,  when 
the  tumors  were  much  smaller. 


Day  64 


Day  72 


Day  81 


Day  29 


Day  43 


Day  52 


Figure  6.  Comparison  of  metastases  detection  ability  of  LET  and  MRI.  Top  row  shows  the  LET 
bioluminescent  signal  overlay  on  a  3D  reconstruction  of  the  animal  surface  obtained  with  external 
illumination.  The  anesthesia  mask  is  seen  over  the  nose  of  the  mouse  in  the  reflected  light  images. 
Bottom  row  shows  representative  slices  from  MRI  scans  covering  the  chest  of  the  mouse,  acquired 
on  a  4.7  T  Varian  scanner. 
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The  figure  shows  four  1mm- thick 
tomographic  slices,  2  for  LET  and  2 
for  MRI. 

LET  shows  three  metastases  on  days 
22  and  43,  two  in  the  left  and  one  in 
the  right  lung.  Using  a  4.5  T  small 
animal  magnet,  MRI  confirms  the 
two  metastases  in  the  left  lung  on  day 
46.  In  the  LET  images  the  maximum 
intensity  of  light  emission  increases 
by  a  factor  of  3-3.5  between  days  22 
and  32. 


Day«  Uay4t> 
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REPORTABLE  OUTCOMES 


In  a  sequential  study  using  the  same  model  the  animal  was  imaged  repeatedly  by  LET  and  MRI 
over  several  weeks  as  the  lung-colonizing  metastases  appeared  and  progressed  in  size  and  number. 
Following  D-luciferin  injection  (450  mg/kg,  SQ)  in  the  anesthetized  mouse,  a  set  of  20  images,  18° 
apart,  was  obtained  using  1  min  exposure  starting  5  min  post-injection.  Data  reconstruction 
provided  a  3-D  model  of  the  lung  tumors.  The  algorithm  was  also  applied  to  images  obtained  using 
light  reflected  from  the  skin  under  external  illumination.  The  imaging  time  varied  from  300 
s/angular  position  for  the  first  imaging  sessions  that  were  characterized  by  low  bioluminescent 
signal,  to  30  s/angular  position  for  the  later  sessions  when  the  bioluminescent  signal  was  large.  MRI 
scans  covering  the  chest  of  the  mouse  were  acquired  on  a  4.7  T  Varian  scanner  using  a  respiratory 
gating  unit.  We  obtained  contiguous  proton  density  weighted  spin-echo  MR  coronal  slices  with  the 
following  parameters:  TE  =  12  ms,  FOV=  3.2  cm  X  6.4  cm,  slice  thickness  =  1  mm,  matrix=  64  X 
128,  4  averages.  The  LET  first  detected  a  bioluminescent  signal  22  days  after  cell  implantation,  and 
revealed  the  growth  and  spread  of  the  lung  metastases  at  weekly  intervals. 

The  tumors  were  first 
detected  by  MRI  46  days 
after  inoculation.  The  smalhl 
tumors  from  the  left  lung 
shown  on  the  LET  images 
were  still  not  observed  by 
MRI.  The  position  and 
relative  sizes  of  the  tumors 
are  consistent  between 
the  two  imaging  modalities. 

However,  LET  was 
able  to  detect  the  tumors 
17  days  earlier 
than  MRI,  when  the  tumors 
were  much  smaller. 


All-: 


01 413-06 
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As  the  tumors  developed,  the  light  intensity  increased  20-fold.  This  would  permit  the  study  of 
smaller  metasases  but  to  study  smaller  foci  and  number  of  cells  it  is  necessary  to  adopt  higher 
magnification  optics.  We  are  investigating  means  to  do  so. 

These  results  validates  LET  against  the  powerful  MRI  3-D  imaging  technique.  The  technique  is 
accurate;  in  addition  it  is  very  sensitive,  and  studies  are  in  progress  to  enhance  this  sensitivity  at 
high  resolution  and  to  validate  the  results  with  increased  statistics. 


MICROSPECT 


Our  project  aims  at  testing  multimodality  imaging,  the  second  modality  being  Single  Photon 
Emission  Tomography.  To  this  end  the  CCDs  are  coupled  to  scintillating  crystals  and  detect  the 
interactions  of  gamma  rays  with  the  crystal.  The  complexity  of  a  practical  device  built  using  this 
technique  are  being  studies  by  a  Ph.D.  student  and  the  results  have  been  submitted  to  the  IEEE  (see 
attachment). 
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25  micron  slit 
1-125 


With  appropriate  optical  means  the  resolution  that  can  be  obtained  is  high.  With  1-125  as  a 
source  and  a  25  pm  slit  the  image  has  a  FWHM  of  60  pm. 

Completion  of  a  small  animal  SPECT  subsystem  for  this  funded  effort  has  been  identified  as  the 
most  important  step  for  the  project.  Although  the  original  proposal  described  imaging  with  iodine- 
125  labelled  agents  to  obtain  the  high  resolution  depicted  above,  the  radiation  burden  to  animals  and 
personnel  necessary  for  obtaining  an  image  is  excessive,  a  decision  was  taken  in  the  fall  of  2006  to 
use  technetium-99m  and  iodine- 123  for  small  animal  imaging.  A  firm  design  goal  in  the  140-159 
keV  range  has  facilitated  several  important  engineering  decisions. 

1.1.1  Selection  of  camera 

Two  in-house  cameras  have  been  developed  for  this  project.  On  paper,  it  was  difficult  to 
determine  if  the  Texas  Instruments  TC253  or  e2v  CCD  67  would  be  a  better  sensor  in  a  lens- 
coupled  gamma  camera.  The  TC253  camera,  developed  using  FPGA  methods  initiated  by  this 
investigator,  completed  in  2006  by  graduate  student  Todd  Soesbe,  and  submitted  for  publication  in 
Fall  2006  by  him,  was  discovered  to  have  an  undocumented  noise  floor  due  to  spurious  charge 
generation  that  was  subsequently  multiplied  in  the  serial  amplification  registers.  A  cooled  e2v 
camera,  developed  by  another  sub-team,  is  not  yet  available  for  testing  and  comparison. 

Rather  that  upgrading  to  the  next  generation  TI  TC  247  sensor,  it  was  decided  to  purchase  4 
commerical  Andor  Luca  cameras  with  this  chip.  This  acquisition  was  not  funded  by  the  project  but 
has  made  it  possible  for  us  to  advance  more  rapidly. 
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Since  multiple  USB  cameras  were  not  supported  by  the  manufacturer’s  software,  we  needed  to 
develop  a  system  for  using  all  cameras  simultaneously.  Software  and  methods  were  developed  for 
eonverting  the  Advanced  Radiological  Sciences  Roeks  4.2  computation  cluster  to  a  Windows  XP 
based  cluster  where  each  node  eontrols  an  emCCD  camera.  The  eomputational  nodes  are  eonnected 
in  a  star  topology  with  a  frontend  through  a  Gigabit  switeh,  and  the  system  ean  be  rapidly 
reeonfigured  from  eomputational  cluster  to  data  acquisition  cluster. 


Gigabit  ethemet  switch 


4  Andor 
Luca 
emCCD 
cameras 


Depending  upon  the  inserted  CDROM,  the  compute  node  will  install  either  Linux  software  for  high 
performance  computing  or  Windows  XP  and  drivers  to  aequire  data  from  the  Andor  Luca.  Images 
are  published  aeross  the  high  speed  network  to  the  frontend  where  proeessing,  analysis,  and  image 
reeonstruction  will  be  developed  in  the  future.  This  high-speed  network  is  private  and  not  aceessible 
from  the  campus  network. 

Selection  of  optics 

Due  to  the  eonservation  of  'etendu  in  opties,  it  is  advantageous  to  minimize  the  field-of-view  on  the 
scintillating  crystal.  Optical  sensitivity  is  hampered  by  demagnification.  Since  the  TC247  sensor  is 
approximately  5  mm  on  a  side,  we  have  estimated  a  2-3  demagnification  factor  will  be  suffieient  to 
obtain  projeetions  from  a  mouse.  Numerous  optieal  sensitivity  experiments  were  performed  with 
various  lens  and  lens  combinations.  We  found  that  a  single  specifie  lens  with  appropriate  rings  to 
modulate  the  field-of-view  was  the  most  sensitive. 

Selection  of  scintillator 

With  a  2-3  cm  field-of-view  finalized,  new  crystal  scintillators  were  acquired.  CsI(Tl) 
remains  the  brightest  scintillator  with  peak  emission  matched  for  solid  state  detectors. 

Selection  of  collimator 

Collimator  design  and  fabrication  remains  a  challenge  for  this  project.  We  have  come  to  believe  that 
pinhole  imaging  with  a  high-intrinsic  resolution  gamma  camera  is  not  the  best  alternative  and  that  it 
is  not  feasible  to  utilize  pinholes  in  the  same  manner  as  they  are  used  with  PMT-based  gamma 
cameras.  Thus  a  parallel-hole  collimator  or  coded  aperture  are  better  suited  for  small  animal 
imaging. 
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CONCLUSIONS 


At  present,  the  four  Luca  Cameras  are  mounted  on  a  test  platform.  Integration  on  an  appropriate 
gantry  -bed  system  would  take  less  than  one  week.  A  new  positioning  gantry  capable  of  supporting 
these  cameras  has  been  designed  and  built. 

We  note  that  if  the  LET  camera  is  used  to  detect  the  skin  of  the  animal  and  the  SPECT  camera 
is  used  simultaneously,  the  perception  of  the  distribution  is  significantly  improved.  In  the  figure 
below,  the  MDDP  distribution  imaged  by  SPECT  identifies  active  regions  of  the  skeleton. 


We  note  that  in  addition  we  can  merge  bioluminescence  with  fluorescence  imaging.  If  a  dextran 
tagged  with  Alexa  Flours  is  injected  in  a  mouse  illuminated  in  white  light,  the  autofluorescence 
presents  a  background  signal  14  times  lower  than  the  signal  produced  by  the  fluor.  In  the  whole 
body  image  the  prominent  light  source  is  the  heart  as  the  polymer  .  Agent  washout  produces  a 
strong  signal  in  the  bladder.is  a  vascular  imaging  agent.  The  were  obtained  in  1  second. 


If  the  same  agent  is  injected  in  a  mouse  in  which  an  MDA-MB231  breast  cancer  tumor  grows  in 
the  mammary  fat  pad,  a  whole  body  image  shows  anatomic  details  as  well  as  the  tumor;  a  higher 
magnification  image  reveals  the  viable  tumor  against  a  darker  halo  around  the  tumor,  produced  by 
absorption  of  the  light  used  to  induce  fluorescence  by  blood  vessels  and  by  a  dense  population  of 
extravasated  red  blood  cells. 
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Images  in  fluorescence  of  a  tumor-bearing  mouse  and  perfusion  of  tumor  compared  to 
contralateral  mammary.  The  dark  ring  surrounding  the  tumor  is  formed  by  necrotic  tissue,  red  blood 
cells,  vascular  tissue. 

Returning  to  SPECT,  we  have  participated  with  our  colleagues  in  evaluating  a  new  imaging 
agent,  Clioquinol.  This  agent  targets  areas  with  high  transition  metal  concentrations.  While  it  was 
originally  developed  to  image  plaque  in  Alzheimer’s  disease  it  is  also  effective  in  breast  cancer. 
Thus  we  show  below  that  the  agent  stains  quite  effectively  a  lymph  node  metastasis  of  breast  cancer 


Lymph  node  metastasis  of  CA 


The  figure  shows  the  radioactivity  distribution  of  '^^I-CQ  in  a  tissue  slice  of  lymph  node 
together  with  a  standard  tumor  stain  in  an  adjacent  slice.  Intense  uptake  is  shown  in  red,  less  intense 
in  blue. 


Storage  phosphor  screen  images  of  BALB/c  nude  mice  injected  with  '^^I-CQ:  tumor-bearing 
mouse  at  right  (MDA-435  Human  BrCa)  vs  control  at  left.  The  control  mouse  shows  a  diffuse 
distribution  but  the  tumor-bearing  mouse  shows  accumulation  in  the  flank. 
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EMCCD  based  Gamma  Camera  for 
Preclinical  SPECT  Imaging 
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Abstract —  The  electron-multiplying  charge-coupled  device 
(EMCCD)  offers  improved  quantum  efficiencies  (40  to  95%)  over 
a  broader  range  of  wavelengths  (400  to  900  nm)  and  a  higher 
intrinsic  resolution  (<  100  /xm  using  photon  counting)  when 
compared  to  photomultiplier  tubes.  The  electron  gain  achieved 
in  the  multiplication  register  of  an  EMCCD  effectively  reduces 
the  readout  noise  to  less  than  1  electron/pixel,  making  them 
sensitive  to  single  photoelectrons.  Our  prototype  camera  uses  the 
Texas  Instruments  Impactron™  EMCCD  model  TC253SPD-B0 
(7.4  /xm  square  pixels)  which  is  cooled  under  vacuum  to  -50  °C 
using  a  four  stage  Peltier  and  liquid  heat  exchanger.  Shuttered 
lens-coupling  is  used  to  image  the  optical  light  from  a  3  mm 
thick  monolithic  CsI(Tl)  crystal.  Precise  clocking  for  the  EMCCD 
is  provided  by  a  National  Instruments  FPGA  controller  (PCI- 
781  IR)  and  Lab  VIEW  EPGA  module  (version  8.0).  A  custom 
built  electronics  box  contains  the  clock  driver  circuitry  and  16-bit 
video  board  for  digital  conversion  of  the  video  signal.  Our  optical 
coupling  method  differs  from  other  EMCCD  SPECT  systems 
by  using  lenses  rather  than  fiber  optic  bundles  for  transfer. 
TC253  characterization  tests  measured  a  maximum  gain  just  over 
1000  X,  dark  current  rate  of  0.14  e/p/s,  read  noise  of  18.2  e/s,  and 
spurious  charge  generation  of  4  e/p.  A  light  integration  intrinsic 
resolution  of  110  /xm  FWHM  was  measured.  Light  integration 
images  of  a  line  phantom  using  a  single  pinhole  collimator  were 
used  for  SPECT  reconstruction.  We  found  the  relative  high 
spurious  charge  generation  and  low  quantum  efficiency  of  the 
TC253  made  it  incapable  of  photon  counting  for  low  energy 
sources  using  lens  coupling. 

Index  Terms — SPECT,  preclinical,  small  animal,  EMCCD, 
gamma  camera,  pinhole  collimator,  high-resolution,  high- 
sensitivity. 

1.  Introduction 

OVER  the  past  decade  the  demand  for  high-resolution 
preclinical  (small  animal)  SPECT  imaging  has  greatly 
increased.  This  increase  is  due  to  the  development  of  new 
radiolabeled  antibodies,  receptor  ligands,  and  other  radiophar¬ 
maceuticals  for  both  imaging  and  therapy  [1].  The  availability 
of  genetically  modified  mice  and  in  vivo  molecular  imaging  for 
studying  gene  expression  and  regulation  in  human  disease  has 
also  increased  the  demand  for  high-resolution  SPECT  [2]-[5]. 
One  drawback  of  previous  small  animal  radiopharmaceutical 
imaging  methods,  such  as  autoradiography,  was  that  several 
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animals  had  to  be  sacrificed  in  order  to  achieve  significant 
statistics  since  each  animal  could  only  be  measured  once  [5]. 
With  the  development  of  high-resolution  preclinical  SPECT, 
an  animal  can  be  measured  several  times  throughout  a  study, 
thereby  allowing  each  animal  to  act  as  its  own  control.  This 
leads  to  fewer  animals  needed  per  study,  which  helps  to  reduce 
the  overall  cost  of  the  experiment. 

The  small  field  of  view  needed  to  image  animals  such  as 
mice  and  rats  (?^  10  cm  x  10  cm)  allows  for  the  use  of  pinhole 
collimation  to  achieve  sub-millimeter  resolution  (<  1  mm 
EWHM).  Yet  the  pinhole  diameter  required  for  this  resolution 
severely  limits  the  sensitivity  of  the  device.  Several  methods 
for  increasing  the  sensitivity  of  preclinical  SPECT  systems 
are  currently  being  pursued.  These  methods  include  using 
multiple  single  pinhole  detectors  [6]-[10],  multiplexing  several 
pinholes  onto  a  single  detector  with  either  minimal  [11],  [12] 
or  extensive  [13],  [14]  image  overlap  on  the  crystal  surface, 
and  even  using  coded  apertures  commonly  found  in  gamma- 
ray  astronomy  [15],  [16].  Existing  clinical  SPECT  systems 
can  be  modified  with  pinhole  collimators  to  achieve  adequate 
resolutions  for  small  animal  imaging  [6],  [13],  [17],  [18],  but 
compact  dedicated  preclinical  SPECT  systems  can  offer  higher 
resolutions  and  are  small  enough  to  fit  on  a  desktop  [7]. 

Another  method  to  improve  sensitivity  is  to  use  a  charge- 
coupled  device  (CCD)  for  the  gamma  camera  photodetector. 
The  application  of  CCDs  to  preclinical  SPECT  imaging  has 
been  previously  shown  by  our  group  [19].  This  research 
has  been  greatly  accelerated  by  the  recent  availability  of  the 
electron-multiplied  CCD  (EMCCD)  which  was  first  techni¬ 
cally  described  in  1992  [20]  and  made  commercially  available 
in  2001  [21],  [22].  EMCCDs  combine  the  positive  aspects  of 
photomultiplier  tubes  (high-gain  leading  to  a  high  signal-to- 
noise  ratio,  photon  counting,  and  energy  discrimination)  with 
the  positive  aspects  of  CCDs  (high  quantum  efficiencies  over 
a  broad  spectral  range,  high  intrinsic  resolution,  and  small 
form  factor)  making  them  well  suited  for  preclinical  SPECT 
applications. 

Our  goal  is  the  development  and  evaluation  of  an  EMCCD 
based  small  animal  SPECT  device  for  murine  model  stud¬ 
ies  of  human  disease.  Previous  applications  of  EMCCDs  to 
preclinical  SPECT  imaging  mainly  use  fiber  optic  bundles  to 
transfer  the  light  from  a  micro-columnar  CsI(Tl)  crystal  to  the 
imaging  area  of  the  EMCCD  [25]-[27],  [42].  This  requires  the 
bundle  to  be  placed  very  close  to,  if  not  touching,  the  EM¬ 
CCD  surface  which  creates  mechanical,  cooling,  and  vacuum 
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Fig.  1.  A  schematic  showing  the  geometry  of  a  gamma  camera  with  a  single 
pinhole  collimator.  When  ^  is  less  than  1,  a  magnified  image  of  the  source 
object  is  projected  onto  the  crystal  surface. 


complications.  Our  method  differs  by  using  monolithic  CsI(Tl) 
crystal  and  lens  coupling  for  light  transfer.  With  lens  coupling, 
physical  contact  with  the  EMCCD  is  eliminated,  although 
the  transfer  efficiency  is  lower  compared  to  fiber  coupling. 
Our  approach  also  differs  by  using  the  very  low  cost  Texas 
Instruments  TC253  EMCCD  in  our  custom  designed  camera. 
This  paper  describes  the  construction  of  the  camera  and 
electronics,  EMCCD  noise  characterization,  gamma  camera 
intrinsic  resolution,  and  SPECT  reconstruction  based  on  light 
integration  and  lens  coupling. 

II.  Materials  and  Methods 
A.  Spatial  Resolution 

“Pinhole  collimators  provide  the  best  combination 
of  sensitivity  and  resolution  for  small  subjects  that 
are  positioned  a  short  distance  from  the  aperture.” 

Hal  O.  Anger,  1967 

As  seen  in  the  above  quote,  the  advantages  of  using  pinhole 
collimators  to  image  small  source  objects  that  are  placed  close 
to  the  aperture  have  been  known  for  some  40  years  [23].  In 
fact,  the  very  first  gamma  camera  used  pinhole  collimation 
to  image  the  uptake  of  ^^^I  (364  keV  gamma  rays)  in  the 
human  thyroid  gland  [24].  The  setup  for  a  single  pinhole 
gamma  camera  is  shown  in  Eigure  1  where  dg  is  the  effective 
diameter  of  the  pinhole  aperture,  b  is  the  distance  from  the 
source  to  the  pinhole,  a  is  the  distance  from  the  pinhole  to  the 
scintillating  crystal,  and  Ri  is  the  intrinsic  spatial  resolution  of 
the  photodetector  and  scintillation  crystal.  The  total  resolution 
of  the  system  is  given  by: 


where  the  first  quadratic  term  under  the  square  root  is  the 
geometric  resolution  (Rg)  and  the  second  quadratic  term  is 
the  effective  intrinsic  resolution  [23].  When  ^  is  less  than  1, 
a  magnified  image  of  the  source  object  is  projected  onto  the 
crystal  surface.  This  makes  the  effective  intrinsic  resolution 
lower  than  the  intrinsic  resolution,  which  helps  to  improve  the 


total  resolution.  Eigure  2(a)  shows  a  plot  of  Rt,  Rg,  and  ^Ri 
versus  ^  for  dg  =  1.0  mm  and  Ri  =  5.0  mm,  which  correspond 
to  the  typical  pinhole  size  used  in  small  animal  imaging  and 
the  typical  intrinsic  resolution  of  a  clinical  SPECT  gamma 
camera  respectively.  Eor  these  conditions,  it  can  be  seen  that 
the  ^Ri  term  quickly  becomes  the  limiting  term  for  the  total 
resolution  as  ^  increases.  As  mentioned  previously,  clinical 
SPECT  gamma  cameras  can  be  used  for  small  animal  imaging 
but  require  large  magnification  factors  to  adequately  image 
mice  and  rats.  Erom  Eigure  2(a)  it  can  be  seen  that  in  order  to 
achieve  a  total  resolution  of  1.5  mm  with  a  clinical  system,  a 
magnification  of  5  x  (^  =  0.2)  must  be  used.  This  leads  to  large 
detector  surface  areas  and  volumes.  In  Eigure  2(b)  dg  is  again 
1.0  mm  but  Ri  has  been  reduced  by  an  order  of  magnitude  to 
0.5  mm.  It  can  be  seen  that  the  geometric  resolution  is  now  the 
limiting  term  for  the  total  resolution,  and  that  a  magnification 
of  only  2x  (^  =  0.5)  is  needed  to  achieve  a  total  resolution  of 
1.5  mm.  This  greatly  reduces  the  detector  surface  area  leading 
to  a  smaller  detector  volume,  which  is  more  suitable  for  small 
animal  imaging. 


B.  Pinhole  Sensitivity 

With  a  very  low  intrinsic  resolution  as  in  Eigure  2(b), 
the  total  resolution  of  the  system  can  be  further  improved 
by  reducing  dg.  Eigure  2(c)  shows  a  resolution  plot  for 
de  =  0.5  mm  and  Ri  =  0.5  mm.  It  can  be  seen  that  Rg 
is  still  the  limiting  term  for  the  total  resolution,  but  now  a 
magnification  of  2x  (^  =  0.5)  gives  a  total  resolution  of 
less  than  1  mm.  Although  having  sub-millimeter  resolution  is 
attractive,  reducing  dg  from  1.0  to  0.5  mm  severely  reduces  the 
sensitivity  (counts/second  per  source  activity)  of  the  system. 
The  on-axis  geometric  efficiency  for  a  pinhole  is  given  by: 


(2) 


and  represents  the  fraction  of  emitted  photons  that  pass 
through  de  from  an  isotropic  point  source  located  a  distance 
b  from  the  pinhole.  Eor  b^  de  the  square  root  term  in  Equa¬ 
tion  2  can  be  binomially  expanded  to  give  the  approximation: 

de^ 


which  is  the  expression  more  commonly  seen  in  the  litera¬ 
ture  [7],  [17],  [23].  When  reducing  de  from  1.0  to  0.5  mm, 
g  goes  from  0.1%  to  0.02%  respectively  (for  b  =  10  mm). 
Therefore  simply  reducing  de  by  a  factor  of  two  reduces  g, 
and  thus  the  sensitivity,  by  a  factor  of  five. 

One  method  to  increase  the  sensitivity  of  pinhole  collima¬ 
tion  is  to  simply  surround  the  source  object  with  multiple 
pinholes.  The  four  main  configurations  for  multiple  pinholes 
are  shown  schematically  in  Eigure  3,  where  the  first  row 
is  the  number  of  pinholes  per  collimator,  the  second  row 
is  the  pinhole  configuration,  and  the  third  row  represents 
the  projected  image  onto  the  crystal  surface  when  imaging 
a  uniform  planar  source.  These  options  include  surrounding 
the  source  object  with  multiple  single  pinhole  detectors  (a), 
multiplexing  several  pinholes  onto  one  photodetector  such 
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d  =  1 .0  mm,  R.  =  5.0  mm 
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d  =  1 .0  mm,  R.  =  0.5  mm 
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Fig.  2.  A  plot  of  Rt,  Rg,  and  ^Ri  versus  ^  for  (a)  de  =  1.0  mm  and  Ri  =  5.0  mm,  (b)  de  =  1.0  mm  and  Ri  =  0.5  mm,  and  (c)  de  =  0.5  mm  and 
Ri  =  0.5  mm.  ^Ri  is  the  limiting  term  for  (a),  while  Rg  is  the  limiting  term  for  (b)  and  (c). 
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Fig.  3.  A  schematic  showing  the  four  main  configurations  for  multiple 
pinhole  collimation.  The  sensitivity  per  collimator  increases  when  moving 
from  left  to  right. 


that  there  is  no  image  overlap  on  the  crystal  surface  (b), 
multiplexing  with  image  overlap  on  the  crystal  surface  (c),  and 
using  a  coded  aperture  which  has  severe  image  overlap  (d). 
Although  the  sensitivity  increases  when  moving  from  left  to 
right  in  Figure  3,  it  is  not  fully  understood  how  the  resolution 
changes  as  the  image  overlap  increases  and  the  tomographic 
reconstruction  becomes  more  complicated.  For  this  project  we 
intend  to  improve  the  sensitivity  by  multiplexing  9  pinholes 
onto  a  single  EMCCD  gamma  camera  in  a  non-overlapping 
manner  and  surrounding  the  source  with  a  ring  of  up  to  12 
cameras  for  a  total  of  over  100  pinholes. 

C.  EMCCD  Sensitivity 

The  sensitivity  of  a  gamma  camera  system  can  be  further 
improved  by  using  an  EMCCD  to  detect  the  optical  photons 
emerging  from  the  scintillation  crystal.  EMCCDs  have  higher 
quantum  efficiencies  (up  to  95%  with  back-illumination)  and 
broader  spectral  responses  (400  to  900  nm)  when  compared 
to  PMTs.  Also,  intrinsic  resolutions  less  than  100  Rm  have 
been  shown  for  EMCCD  based  gamma  cameras  [25],  [26].  The 
most  important  characteristic  of  an  EMCCD  is  the  improved 
low  light  level  signal-to-noise  ratio  (SNR)  when  compared 
to  conventional  CCDs.  This  is  achieved  by  the  addition  of 
a  second  serial  register,  known  as  the  gain  or  multiplication 
register,  to  the  CCD  architecture  as  shown  in  Figure  4. 


Electrons  are  first  captured  within  the  pixel  potential  wells  of 
the  image  area,  then  the  entire  frame  (for  example  680x500 
pixels)  is  parallel-transferred  to  the  storage  area.  A  single 
pixel  row  is  then  parallel-transferred  from  the  storage  area 
to  the  serial  register.  The  pixel  row  in  the  serial  register  is 
then  serial-transferred  through  the  multiplication  register  to 
the  charge-to- voltage  output  amplifier  of  the  EMCCD.  This 
is  done  for  each  row  until  the  entire  storage  area  has  been 
read.  The  multiplication  register  differs  from  the  normal  serial 
register  by  the  addition  of  a  third  gate  known  as  the  charge 
multiplication  gate  (CMC)  [28].  By  making  the  voltage  of 
the  CMC  very  high  (up  to  15  VDC)  the  electrons  undergo 
impact  ionization  as  they  move  through  the  multiplication 
register,  which  amplifies  the  electric  signal.  Although  the  gain 
per  pixel  step  is  very  low  (1.01  to  1.018)  a  multiplication 
register  of  400  pixels  in  length  can  allow  gains  of  over 
lOOOx  to  be  reached  [29].  Multiplication  of  the  signal  in  the 
charge  domain  before  charge-to-voltage  conversion  effectively 
reduces  the  readout  noise  of  the  output  amplifier  to  less 
than  1  electron/pixel.  This  eliminates  the  noise  fioor  found 
in  conventional  CCDs  making  EMCCDs  sensitive  to  single 
electrons. 

The  total  signal  noise  of  an  EMCCD  is  given  by: 

at  =  v/G2F2(P  +  D  +  5)+a2  (4) 


where  G  is  the  gain,  F  is  the  excess  noise  (Fano)  factor  due 
to  the  gain  process  in  the  multiplication  register,  and  P,  D, 
S,  and  ar  are  the  mean  signals  (electron/pixel)  due  to  the 
photon  (shot),  dark  current,  spurious  charge,  and  read  noises 
respectively  [29],  [30].  The  dark  current  noise  (D)  and  the 
spurious  charge  noise  (S)  can  be  reduced  to  negligible  levels 
with  cooling  and  proper  clocking  respectively  [31].  Also,  for  a 
gain  greater  than  10  x,  F  approaches  ^/2  [32].  Under  these 
conditions  the  SNR  for  a  EMCCD  reduces  to: 


GP 

SNRemccd  = - —  \  F 

at  V  2 

when  G  ^  ar.  This  can  be  compared  to  the  SNR  of  a 
conventional  CCD  under  similar  condition  given  by: 


SNR  CCD 
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Fig.  5.  A  plot  showing  the  EMCCD  and  CCD  SNR  as  a  function  of 
electrons/pixel.  The  EMCCD  SNR  is  much  higher  for  low  light  level  imaging. 


multiplication  register 
400  pixels 


output  amplifier 


Fig.  4.  A  schematic  showing  the  architecture  of  an  EMCCD.  The  electrons 
are  multiplied  by  impact  ionization  as  they  move  through  the  multiplication 
register  before  charge-to-voltage  conversion  in  the  output  amplifier. 

A  plot  of  SNRemccd  and  SNRccd  versus  P  is  shown  in 
Figure  5  for  a  typical  read  noise  =  25  electrons/pixel). 
It  can  be  seen  that  while  the  CCD  has  a  SNR  of  1  at 
25  electrons/pixel,  the  EMCCD  has  a  SNR  that  is  3.5  times 
higher. 


Eig.  6.  A  schematic  showing  the  setup  for  the  prototype  EMCCD  camera. 


D.  Camera  Prototype 

A  schematic  of  the  prototype  system  is  shown  in  Figure  6. 
The  system  consists  of  a  PC  with  a  FPGA  card,  an  external 
electronics  box,  and  the  EMCCD  camera  assembly.  The  Na¬ 
tional  Instruments  PCI-781  IR  (Xilinx  FPGA,  1  million  gates, 
40  MHz  digital  rate,  25  ns  resolution)  produces  the  clock 
pulse  timing  for  the  EMCCD.  These  pulses  are  sent  to  the 
driver  board  inside  the  external  electronics  box  which  converts 
the  TTE  pulses  to  the  proper  voltages  needed  to  drive  the 
EMCCD.  The  driving  pulses  are  then  sent  through  an  R/C 
circuit  located  inside  the  camera  housing  for  shaping  before 
entering  the  EMCCD.  The  video  signal  from  the  EMCCD  is 
first  sent  through  a  preamp  and  then  through  a  video  board. 
The  video  board  uses  correlated  double  sampling  (CDS)  of  the 
video  signal  to  help  reduce  noise  before  digital  conversion  in 
a  16-bit  ADC  (Analog  Devices  AD976A,  200  kHz  max  rate). 
The  digitized  video  signal  is  then  sent  back  to  the  FPGA  which 
uses  a  DMA  FIFO  across  the  PCI  bus  to  store  the  data  on  the 
host  computer  which  is  running  Lab  VIEW  8.0. 

For  the  prototype  camera  we  used  the  Texas  Instruments 
Impactron™  TC253  EMCCD  chip  (TC253SPD-B0).  This 
front-illuminated  chip  (35%  quantum  efficiency  at  560  nm) 
has  7.4  pm  square  pixels  (100%  fill  factor)  in  a  680x500, 
1/3  inch  format.  The  TC253  has  a  nominal  read  noise  of 


29  electrons/pixel  RMS  at  video  rate  readout  (12.5  MHz 
clock  frequency).  It  also  has  low  dark  current  generation 
0.1  electron/pixel/sec  at  -40  °C)  and  a  built  in  single  stage 
Peltier  element  capable  of  cooling  the  chip  down  to  -10  °C. 
Figure  7(a)  shows  a  close  up  photo  of  the  TC253.  The  TC253 
is  mounted  onto  an  Aluminum  cold-finger  which  sits  on  top  of 
a  four- stage  Peltier  element.  Low  thermal  conductivity  wires 
connect  the  TC253  to  the  R/C  and  preamp  circuitry  with  all 
camera  electronics  being  contained  in  the  vacuum  housing. 
Figure  7(b)  shows  an  exploded  view  of  the  camera  assembly 
showing  the  lens,  shutter,  vacuum  window,  EMCCD  housing, 
and  a  choice  of  either  air  or  liquid  heat  exchanger.  A  picture 
of  the  assembled  camera  is  shown  in  Figure  7(c). 

HI.  Results 

High  quantum  efficiency  back-illuminated  EMCCDs  such 
as  the  e2v  CCD97  currently  cost  between  6k  to  7k  USD.  Off- 
the-shelf  camera  systems  containing  these  EMCCDs  such  as 
the  PhotonMAX  512B  cost  around  30k  USD.  To  increase  the 
sensitivity  of  our  preclinical  SPECT  system,  we  would  like  to 
surround  the  animal  with  a  ring  of  multiple  cameras.  There¬ 
fore,  in  order  to  reduce  the  cost  per  camera  to  a  reasonable 
level,  we  have  opted  to  design  and  build  our  own.  By  using 
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Fig.  7.  (a)  A  photo  showing  the  front  face  of  the  Texas  Instruments  Impactron™  TC253  EMCCD  mounted  to  the  Aluminum  cold  finger.  The  image  area 

(680  X  500  pixels)  measures  5.0  mm  x  3.7  mm.  (b)  An  exploded  photo  of  the  camera  assembly,  (c)  The  assembled  camera  with  air  heat  exchanger. 


the  low-cost  front-illuminated  TC253  EMCCD  (500  USD), 
the  PCI-781  IR  FPGA  card  with  Lab  VIEW  module  (2k  USD), 
and  custom  built  electronics  (500  USD),  the  cost  per  camera 
is  reduced  by  an  order  of  magnitude  to  3k  USD.  Details  of 
the  camera  system,  followed  by  TC253  noise  characterization, 
and  initial  light  integration  imaging  results  are  presented  next. 
Although  the  system  has  a  gain  of  over  lOOOx,  low  read  and 
dark  noises,  and  good  intrinsic  resolution,  the  relative  high 
spurious  charge  generation  and  low  quantum  efficiency  of  the 
TC253  make  it  incapable  of  photon  counting  for  both 
and  ^^^I  using  monolithic  CsI(Tl)  crystals  with  lens  coupling. 

A.  EMCCD  Characterization 

Upon  completing  the  prototype  camera,  we  first  character¬ 
ized  the  EMCCD  by  measuring  the  gain,  dark  current  accu¬ 
mulation  rate,  total  read  noise,  and  spurious  charge  generation 
to  verify  that  they  were  within  the  specified  operating  range 
of  the  TC253.  Methods  and  results  for  the  measurements  are 
described  here.  All  measurements  were  taken  at  -40  °C. 

1)  Gain:  The  gain  versus  positive  charge  multiplication 
gate  voltage  (-fCMG)  of  the  TC253  was  measured  using 
lens  coupling  to  image  a  controlled  light  source.  The  light 
source  consisted  of  a  3  mm  green  LED  embedded  in  a 
10  mm  X  15  mm  Tefion  cylinder  to  create  a  small  diffuse 
uniform  source  that  was  viewed  on  end.  With  the  EMCCD 
gain  set  to  1  (-fCMG  =  7.0  V),  the  current  to  the  LED  was 
varied  until  a  SNR  of  about  5  was  achieved  for  an  average 
of  1000,  33  msec  exposures.  The  difference  between  the  dark 
reference  pixels  and  peak  signal  was  measured  using  a  line 
profile  of  the  averaged  image.  This  measurement  was  repeated 
as  the  -fCMG  voltage  was  increased,  with  the  increase  between 
the  dark  and  peak  pixels  being  taken  as  the  gain.  A  plot  of 
gain  versus  -i-CMG  is  shown  in  Figure  8.  It  can  be  seen  that 
the  gain  starts  to  increase  past  11  V  and  rises  to  over  1000  at 
14  V. 

2)  Dark  Current:  Dark  current  arises  from  thermal  gener¬ 
ation  of  electrons  within  the  silicon  lattice  of  the  EMCCD. 
These  electrons  are  captured  by  the  pixel  potential  wells  and 
contribute  to  the  non-photon  related  signal  noise  (see  Equa¬ 
tion  4).  Dark  current  generation  is  proportional  to  exposure 
time,  independent  of  the  amount  of  light  incident  on  the  EM¬ 
CCD,  and  can  be  significantly  reduced  by  cooling  [33].  The 
dark  current  accumulation  rate  for  the  TC253  was  measured 
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Fig.  8.  A  plot  of  EMCCD  gain  versus  positive  charge  multiplication  voltage 
(+CMG)  for  the  TC253.  A  gain  of  over  lOOOx  is  reached  at  14  V. 

by  taking  N  exposures  at  a  fixed  dark  exposure  time.  The 
mean  and  standard  deviation  of  the  N  exposures  were  then 
calculated  for  each  pixel.  The  ensemble  average  of  the  pixel 
means  over  the  active  EMCCD  area  was  taken  as  the  dark 
current  count  (electrons/pixel),  while  the  ensemble  average 
of  the  pixel  standard  deviations  was  used  to  calculate  the 
error  (±-^).  This  measurement  was  then  repeated  for  several 
different  dark  exposure  times.  For  these  measurements,  the 
gain  was  set  to  1  and  the  EMCCD  temperature  was  -40  °C. 
A  plot  of  exposure  time  versus  dark  current  count  is  shown 
in  Figure  9.  The  fitted  slope  gives  0.14  electrons/pixel/second 
which  agrees  with  the  dark  current  rate  given  in  the  TC253 
datasheet. 

3)  Read  Noise:  The  main  source  of  read  noise  in  the  system 
arises  from  the  charge-to-voltage  conversion  in  the  EMCCD 
output  amplifier.  Assuming  that  the  spurious  charge  has  been 
reduced  to  a  negligible  level,  then  for  a  zero  second  dark 
exposure  (P  =  D  =  0)  at  a  gain  of  1,  Equation  4  reduces 
to: 

(Jt  =  =  (Jr  (7) 

SO  that  the  total  noise  is  equal  to  the  total  read  noise.  The  read 
noise  of  the  system  was  measured  by  taking  1000  zero  second 
exposures  at  a  gain  of  1  and  calculating  the  standard  deviation 
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Fig.  9.  A  plot  of  dark  current  count  (electrons/pixel)  versus  dark  exposure 
time.  The  fitted  slope  gives  the  dark  current  accumulation  rate  for  the  EMCCD 
at  -40  °C. 


for  each  pixel.  The  ensemble  average  of  the  pixel  standard 
deviations  over  the  active  EMCCD  area  was  then  multiplied 
by  the  system  video  gain  (0.57  electrons/ADU)  to  give  a  read 
noise  of  18.2  electrons/pixel  RMS.  The  measured  read  noise 
was  lower  than  the  read  noise  stated  in  the  TC253  datasheet 
(29  electrons/pixel  RMS)  due  to  the  200  kHz  maximum 
readout  rate  of  the  current  ADC.  Although  fine  for  initial 
testing,  a  faster  ADC  will  be  needed  to  achieve  video  rate 
readout  (12.5  MHz).  This  will  increase  the  readout  noise  which 
is  proportional  to  the  readout  speed. 

4)  Spurious  Charge:  Along  with  dark  current,  non-photon 
related  electrons  known  as  spurious  charge  (or  clock  induced 
charge)  can  be  generated  during  the  readout  of  the  EMCCD. 
Spurious  charge  arises  from  the  rapid  inversion  of  the  clocking 
potentials  and  tends  to  increase  with  faster  readout  speeds  and 
lower  temperatures  [33].  Methods  to  reduce  spurious  charge 
include  using  a  slow  parallel-transfer  from  the  image  area  to 
the  storage  area,  setting  the  rise  and  fall  times  of  the  clocking 
pulses  to  about  200  ns,  and  using  the  minimum  clocking 
potentials  necessary  for  efficient  charge  transfer  [26],  [31]. 
Of  these  three  methods,  the  greatest  observed  reduction  in 
spurious  charge  for  the  TC253  was  seen  by  minimizing  the 
clocking  potentials  [34]. 

Spurious  charge  generation  as  a  function  of  -fCMG  was 
measured  for  the  TC253  by  taking  100  zero  second  dark 
exposures  for  a  fixed  -fCMG  voltage.  The  ensemble  average  of 
the  standard  deviations  for  each  pixel  was  then  calculated  and 
taken  as  the  total  noise  (at).  Eor  a  zero  second  dark  exposure 
(P  =  D  =  0)  Equation  4  can  be  rewritten  as: 

*  -  p2Q2  (8) 

By  using  at  from  above,  a^  and  G  from  the  previous  read 
noise  and  gain  measurements,  and  knowing  that  F  \f2  at 
high  gain,  the  spurious  charge  generation  was  calculated  to 
be  3  to  4  electrons/pixel  at  gains  from  500  x  to  lOOOx.  This 
spurious  charge  generation  could  not  be  further  reduced  by 
the  minimizing  methods  listed  above  and  appears  to  be  an 
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Fig.  10.  A  schematic  of  the  EMCCD  SPECT  setup.  Two  lenses  are  used 
to  improve  the  optical  transfer  of  the  scintillation  light  to  the  EMCCD.  The 
lens  focal  plane  can  be  moved  through  the  crystal. 


inherent  property  of  the  TC253,  largely  caused  by  the  use  of 
virtual-phase  electrode  technology  [22].  This  spurious  charge 
generation  is  therefore  the  true  noise  fioor  for  the  TC253 
chip  making  it  incapable  of  single  photoelectron  sensitivity. 
Consequences  of  this  measurement  are  further  discussed  in 
Section  IV.  As  a  comparison,  this  same  measurement  was  done 
with  a  high-quality  commercially  available  EMCCD  camera 
(Princeton  Instruments/ Acton  PhotonMAX  512B  with  the  e2v 
CCD97  EMCCD)  which  yielded  a  spurious  charge  generation 
of  less  than  0.01  electrons/pixel. 

B.  Intrinsic  Resolution 

A  schematic  of  the  lens-coupled  SPECT  setup  is  shown 
in  Eigure  10.  Two  1-inch  format  lenses  (Schneider  Xenon 
0.95/25)  were  coupled  front-to-front  in  order  to  improve  the 
optical  transfer  to  the  EMCCD.  Both  lenses  had  an  f-stop 
of  0.95  with  the  front  lens  focus  set  to  infinity  and  the  rear 
lens  focus  being  adjustable.  A  3  mm  thick  monolithic  CsI(Tl) 
crystal  was  then  placed  5  to  10  mm  in  front  of  the  lenses  so 
that  the  shallow  focal  plane  (less  than  1  mm  deep)  could  be 
moved  through  the  crystal.  The  two  coupled  lenses  image  a 
field  of  view  (EOV)  on  the  crystal  that  is  approximately  the 
same  size  as  the  image  area  of  the  EMCCD  (?^  5  mm  x  5  mm). 
When  compared  to  a  single  lens  imaging  the  same  EOV, 
we  found  the  optical  transfer  when  using  two  lenses  to  be 
improved  by  a  factor  of  two. 

The  intrinsic  resolution  (Ri)  for  the  above  setup  was  mea¬ 
sured  using  a  30  pm  slit  Tungsten  collimator  and  a  ^^’^Tc 
source.  The  activity  viewed  by  the  slit  was  approximately 
50  pC\  (1.9  MBq).  The  line  spread  function  for  a  5  minute 
exposure  with  the  gain  set  to  200  x  was  acquired.  The  acquired 
image  and  1  mm  wide  line  profile  (taken  perpendicular  to  the 
slit)  are  shown  in  Eigure  11.  Eitting  the  line  spread  function 
profile  with  a  Gaussian  distribution  gave  a  EWHM  of  110  pm 
which  was  taken  as  the  intrinsic  resolution. 

C.  SPECT  Images 

Eor  SPECT  imaging,  the  Tungsten  slit  collimator  was  re¬ 
placed  with  a  single  pinhole  Lead  collimator.  The  pinhole 
diameter  was  0.5  mm  with  an  aperture  angle  of  90  degrees. 
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Fig.  II.  The  image  and  line  profile  produced  by  a  30  fim  slit  collimator 
and  99"^  Tc  source.  The  exposure  time  was  5  minutes  with  a  gain  of  200  x. 


A  small  glass  capillary  tube  was  filled  with  aqueous 
solution,  mounted  to  a  step  motor,  and  placed  in  front  of  the 
pinhole.  The  3  cm  long  capillary  had  an  outer  diameter  of 
1.5  mm,  an  inner  diameter  of  0.8  mm,  and  a  linear  activity  of 
10  /iCi/mm  (0.37  MBq/mm).  Projection  data  were  then  taken 
as  the  phantom  was  rotated  (20  projections,  18°  step)  with  a 
3  minute  exposure  time  for  each  projection.  The  EMCCD  gain 
was  set  to  1000 x  at  a  temperature  of  -50  °C.  Figure  12  shows 
a  sample  projection  of  the  line  source,  with  the  subject  FOV 
measuring  6.0  mm  x  4.5  mm. 

Image  reconstruction  was  based  on  the  Maximum  Fikeli- 
hood  method  (MFEM)  [35]  which  was  improved  with  the 
listmode  high  resolution  algorithm  for  system  matrix  modeling 
(MESM)  [36],  [37].  A  convolution  procedure  based  on  the 
discrete  Fast  Fourier  Transform  (FFT)  was  used  for  resolution 
modeling.  If  photon  scattering  is  neglected,  the  blurring  com¬ 
ponent  of  the  matrices  which  account  for  the  finite  resolution 
effects  can  be  represented  as  a  set  of  shift-invariant  kernels 
(details  given  in  [37],  [38]).  Figure  13(a)  shows  a  transverse 
slice  of  the  MFSM  reconstruction  for  the  line  phantom.  The 
total  FOV  (20  mm^)  was  subdivided  into  0.25  mm^  voxels. 
Approximately  300,000  events  were  reconstructed  into  the 
voxel  array.  A  Gaussian  function  (0.4  mm  kernel)  was  used 
to  model  the  resolution  during  the  10  iterations.  The  FWHM 
of  the  Gaussian  fitted  line  profile  in  Figure  13(b)  measures 
0.96  mm.  The  slight  increase  in  the  reconstructed  diameter  is 
due  to  the  total  resolution  of  the  setup  being  limited  by  the 
geometric  resolution  {Rg)  which  was  1  mm. 

IV.  Conclusion  and  Discussion 
A.  Summary 

A  cost  effective  prototype  EMCCD  based  gamma  camera 
has  been  developed  and  evaluated  for  use  in  preclinical  SPECT 
imaging.  A  four  stage  Peltier  and  liquid  heat  exchanger  was 
used  to  cool  the  EMCCD  down  to  -50  °C.  Dual  f/0.95 


Fig.  12.  A  pinhole  projection  image  of  a  filled  line  phantom.  The 

subject  FOV  measures  6.0  mm  x  4.5  mm  with  an  exposure  time  of  3  minutes. 


(a)  Z,  mm  (b)  Z,  mm 

Fig.  13.  SPECT  reconstruction  of  the  line  source  (a)  with  Gaussian  fitted 
line  profile  (b). 


lens  coupling  was  used  to  image  the  scintillation  light  from 
a  3  mm  thick  monolithic  CsI(Tl)  crystal.  A  gain  of  over 
1000 X  was  measured  for  the  TC253  EMCCD.  The  dark 
current  accumulation  rate  and  read  noise  were  within  the 
specified  operating  range  for  a  temperature  of  -40  °C  and 
readout  rate  of  200  kHz.  A  spurious  charge  generation  of  3 
to  4  electrons/pixel  was  measured  at  high  gain  values.  This 
intrinsic  spurious  charge  could  not  be  further  reduced  using 
established  minimizing  techniques  and  represents  the  noise 
floor  for  the  TC253.  An  intrinsic  resolution  of  1 10  jim  FWHM 
was  measured  using  a  30  jim  slit  collimator  and  source. 

Single  pinhole  collimation  was  used  to  take  SPECT  images  of 
a  small  ^^"^Tc  line  phantom.  The  MFEM  reconstructed  source 
size  was  slightly  larger  than  the  true  source  size  due  to  the 
geometric  resolution  of  the  setup. 

B.  Conclusions 

Although  the  prototype  EMCCD  gamma  camera  was  capa¬ 
ble  of  SPECT  imaging  using  light  integration,  it  was  unable  to 
perform  photon  counting  for  ^nd  ^^^I  with  the  current 

optical  transfer  method.  This  can  be  explained  by  estimating 
the  number  of  electrons  produced  in  the  TC253  per  CsI(Tl) 
scintillation  event.  The  luminosity  of  CsI(Tl)  is  54  pho- 
tons/keV  with  a  peak  wavelength  of  560  nm.  Approximately 
10%  of  these  scintillation  photons  will  escape  the  crystal  via 
the  non-internally  refiected  light  cone  [39].  Upper  limits  for 
dual  lens  optical  coupling  efficiency  (with  an  object  to  image 
ratio  of  1:1)  have  been  estimated  to  be  5%  [40],  [41], 
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and  the  quantum  efficiency  of  the  TC253  is  35%  at  560  nm. 
This  gives  (140x54x0.1x0.05x0.35)=  13  electrons/event  for 
99mTc  and  (27x54x0.1x0.05x0.35)  =  3  electrons/event  for 
When  these  signals  are  spread  over  a  3  x  3  pixel  area  they 
are  easily  lost  within  the  3  to  4  electrons/pixel  spurious  charge 
noise  of  the  TC253.  Single  photoelectric  events  were  observed 
using  a  higher  energy  ^^^Cs  (662  keV)  source  where  the 
increased  number  of  scintillation  photons  results  is  more  than 
3  to  4  electrons/pixel/event.  Also,  when  the  same  scintillation 
crystal  and  lens  setup  were  used  with  the  PhotonMAX  512B 
camera  (e2v  CCD97  EMCCD,  95%  quantum  efficiency  at 
560  nm)  single  photoelectric  events  were  easily  observed  for 
both  99"^Tc  and  sources.  These  results  are  to  be  presented 
in  a  subsequent  publication. 

C.  Further  Work 

One  possible  solution  to  allow  photon  counting  with  the 
TC253  is  to  improve  the  optical  transfer  of  the  scintillation 
photons  to  the  EMCCD  by  using  a  micro-columnar  CsI(Tl) 
crystal  with  a  fiber-optic  bundle  instead  of  monolithic  crystal 
with  lens  coupling.  Eiber-optic  coupling  offers  an  order  of 
magnitude  improvement  in  the  coupling  efficiency  (>  50% 
with  a  1:1  image  to  object  ratio)  compared  to  lens  cou¬ 
pling  [39].  This  method  of  photon  counting  has  been  previ¬ 
ously  established  [25],  [26],  [42]  and  would  allow  for  energy 
windowing,  a  higher  intrinsic  resolution,  and  a  larger  crystal 
EOV  to  be  observed.  Another  solution  could  be  to  change 
from  the  Texas  Instruments  TC253  to  the  TC247  EMCCD.  The 
TC247  has  a  higher  quantum  efficiency  than  the  TC253  (50% 
versus  35%  at  560  nm)  and  does  not  appear  to  be  noise  limited 
by  inherent  spurious  charge  generation  [43].  Although  the 
TC247  has  a  larger  pixel  size  than  the  TC253  (10  jim  x  10  jim 
versus  7.4  jim  x  1 A  yum)  it  does  not  have  a  100%  fill  factor 
which  could  effect  the  sensitivity.  These  two  options  will  be 
explored  with  the  goal  of  making  our  EMCCD  gamma  camera 
capable  of  both  photon  counting  and  energy  discrimination  for 
99mTc  and  sources. 

D.  Discussion 

EMCCD  images  formed  by  counting  individual  photo¬ 
electric  events  within  the  scintillation  crystal  have  improved 
intrinsic  resolutions  when  compared  to  light  integration  im¬ 
ages  [25].  Although  resolutions  of  less  than  100  /im  have 
been  measured,  it  is  not  yet  understood  how  to  take  full 
advantage  of  this  metric  in  preclinical  SPECT  imaging.  Eor 
example,  when  using  a  0.5  mm  pinhole  collimator  with  an 
intrinsic  resolution  of  100  /im,  the  total  spatial  resolution 
is  still  limited  by  the  geometric  resolution  of  the  pinhole 
(see  Equation  1).  The  total  resolution  can  be  improved  by 
reducing  the  pinhole  diameter,  but  this  has  negative  effects 
on  the  sensitivity  as  previously  shown.  Eigure  14  shows  a 
resolution  plot  for  7?^  =  dg  =  100  /im.  It  can  been  seen 
that  for  ^  =  1  (no  image  magnification)  a  total  resolution  of 
about  225  /im  is  achieved.  It  can  also  be  seen  that  a  total 
resolution  of  less  than  1  mm  is  still  reached  with  -  =  5, 
where  the  projected  image  on  the  crystal  surface  is  smaller 
than  the  subject  object.  This  means  that  for  a  subject  EOV 


d  =  0.1  mm,  R.  =  0.1  mm 
e  I 


Fig.  14.  A  plot  of  Rt,  Rg,  and  ^Ri  versus  ^  for  de  =  0.1  mm  and 
Ri  =0.1  mm.  It  can  be  seen  that  a  total  resolution  of  less  than  1.0  mm  is 
still  achieved  at  -  =5. 

a 

of  10  mm  X  10  mm,  only  2  mm  x  2  mm  is  needed  on 
the  crystal  surface.  Therefore,  one  hundred  0.1  mm  diameter 
pinholes  could  be  multiplexed  onto  a  20  mm  x  20  mm 
collimator  with  no  image  overlap  on  the  crystal  surface.  Yet  the 
combined  sensitivity  of  these  one  hundred  pinholes  would  be 
equivalent  to  a  single  1.0  mm  diameter  pinhole.  It  is  possible 
that  micro-channel  parallel  hole  [44],  [45]  or  multiple  slit  [46] 
collimators  might  better  use  the  high  intrinsic  resolution  of  an 
EMCCD  based  gamma  camera  to  find  a  good  balance  between 
sensitivity  and  resolution  in  preclinical  SPECT  imaging. 
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